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Abstract 



Lepton flavour non-conserving processes are examined in the context of uni- 
fied models with {7(l)-family symmetries which reproduce successfully the 
low-energy hierarchy of the fermion mass spectrum and the Kobayashi - 
Maskawa mixing. These models usually imply mixing effects in the super- 
symmetric scalar sector. We construct the fermion and scalar mass matrices 
in two viable models, and calculate the mixing effects on the \x — > ej, \x — > 3e 
and t — > /i7 rare decays. The relevant constraints on the sparticle mass spec- 
trum as well as the role of various MSSM parameters are discussed. 
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One of the most dramatic consequences of supersymmetric extensions 
of the Standard Model (SM) is the appearance of new sources of flavour 
violations jl], @, 0]. Supersymmetric partners of gauge, fermion and scalar 
fields generate new types of flavour violating diagrams at the one-loop level, 
which enhance considerably the various rare processes. 

Flavour non- conserving processes may still be relatively suppressed if the 
matrices of the supersymmetric partners of fermions, i.e. those of scalar 
quarks and scalar leptons, are diagonal in flavour space. It is widely be- 
lieved however, that a realistic spectrum for the fermion mass matrices can 
be obtained when additional symmetries discriminate between the various 
families of the known fermion fields. Such symmetries imply also a non- 
trivial structure for the corresponding scalar mass matrices. Rare processes, 
being sensitive to these changes, usually lead to hard violation of flavour. 

In this work, we compute the branching ratios for lepton flavour violating 
decays in realistic models whose fermion and scalar mass textures are ob- 
tained by £7(l)-family (£7(1)/-) symmetries. In our analysis we choose both, 
symmetric and non-symmetric fermion textures by appropriate selection of 
the U(l)f fermion charges. The choice of the lepton sector to test the pre- 
dictability, and possibly the viability, of U(l) /-models is ideal as there are 
no large uncertainties (unlike the case of the quark sector where large am- 
biguities enter due to poor knowledge of hadronic matrix elements) in the 
calculations. In a previous work 0], we have given some estimates for the 
/i — > e7 decay is the context of a simple t/(l)/-model in the case of small 
tan/3 regime. Here, we extend our previous analysis and examine cases for 
large and small tan (3 and various values of the gravitino and gaugino masses. 
We use two-loop analysis for gauge couplings and take into account thresh- 
old effects to calculate the sparticle spectrum used to construct the relevant 
scalar and fermion mass matrices entering the above processes. We find that 
the non-observation of lepton flavour violating processes put rather strong 
lower limits on the sparticle mass spectrum, in particular when tan /3 is large. 

A wide class of models, which naturally incorporate flavour non-diagonal 
scalar mass matrices, arises in string scenarios where the usual gauge sym- 
metry is accompanied by a number of U(l) factors, the latter playing the 
role of family symmetries. The fermion mass textures of the above models 
are dictated by the particular charges of the particles under the U{l)f sym- 
metries, the specific flat direction which has been chosen as well as the string 
selection rules and other string symmetries. In general, there are only few 
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tree-level couplings in the superpotential (usually only those responsible for 
the top, bottom and tau masses), while all other fermion mass entries are 
supposed to be generated by higher non-renormalizable terms. 

Once the flat directions and the U(l)f charges of a particular model have 
been fixed, the scalar mass matrix structure may also be easily computed 
through the Kahler function 

£ = /c + io g |nf (i) 

where W is the superpotential and /C has the general form 

K = - \og{S + S) - Y, K log(T„ + f n ) + Z i3 (T n , f n )Q t Q 3 + ■■■ (2) 

with Qi being the matter fields, S the dilaton, whereas T n are the other 
moduli fields. The scalar mass matrices are determined by Zq and W. The 
form of the Zq function is dictated by the modular symmetries and depends 
on the moduli and the modular weights of the fields. Thus, at the tree level, 
the diagonal terms are the only non-zero entries in the scalar mass matrices. 
Higher order terms allowed by the symmetries of the specific model fill in the 
non-diagonal entries. In what follows, we will explore the flavour violating 
processes in two different models which give realistic fermion mass spectrum. 
In particular, we will calculate the fj, — ■> ej, r — > /ry and ji — > 3e processes in 
the context of supersymmetric models whose low energy theory is the MSSM 
model augmented by a U(l) family symmetry. One of them is using a charge 
assignment where symmetric mass matrices appear, while the other assumes 
[/(l)-charges which lead to non-symmetric textures. 

We start with some preliminary remarks about the sources of flavour vi- 
olations and set our notation and conventions. After the breaking of some 
possible unified symmetry to that of the Standard Model, the Yukawa inter- 
action, which appear in the superpotential and violate lepton flavour, is 

W = e cT X e £H + ■■■ (3) 

where £ is the left lepton doublet, e is the right singlet lepton, H is the higgs 
doublet and A e represents the Yukawa coupling matrix in flavour space. In 
addition, soft supersymmetry breaking terms generate mass matrices for the 
charged slepton fields, denoted by rhe, fh eR . The Yukawa and soft scalar 
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mass-squared matrices are diagonalized by unitary transformations 

A e = V* R \ s e Vl (4) 
U e (mf) 5 Uj, and m\ R = U eR (mi R ) * U\ R (5) 



m 



where 5 denotes diagonal. The lepton mass eigenstates (eig) are related to 
the weak eigenstates (w) by 



-eig 



Vle c w (6) 



The charged slepton mass-squared matrix is a 6 x 6 matrix, built up from 
the the two 3x3 left mf and right ra^ soft ones, as well as the off-diagonal 
submatrix which has the form 



( m LR) = m i{^ tan /5 + A t) (7) 

where mi = X e v cos/3/a/2 is the charged lepton mass matrix, Ai and /i are 
the trilinear and higgs mixing parameters in the superpotential, tan/? the 
ratio of the two higgs vev's and v = 246GeV. 

To calculate the mixing effects in the amplitudes, we work in the basis 
where the fermion mass matrix is diagonal, 

(e cT ) w md w = (e cT ) w V^m s e Vle w = (e c ) T eig m 5 e e eig (8) 

In this basis, the off-diagonal term (^) is written 

m^Otan/3 + A e ) = V^mf (// tan (3 + A t )Vl (9) 

This defines uniquely the scalars in the basis where the fermions are in their 
mass-eigenstates 



e* 



R = (vUe R )) w = V£(e R r w (10) 

i' = viL. (ii) 

The soft terms for right and left charged sleptons must be written in the 
same basis. The 6x6 matrix then takes the form 



VimjV L ((A e + fi tan 0)mf )* 



(12) 
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A well known result in the context of the non-supersymmetric standard 
model is the conservation of lepton flavour in the case of zero neutrino masses, 
while in the case of massive, non-degenerate neutrinos, the amount of lepton 
flavour violation is proportional to the factor Am^/M^, ||, which highly 
suppresses all relevant processes. When supersymmetry enters the game, 
the whole scene changes completely. Even in the absence of right handed 
neutrinos, flavour violations could occur via the exchange of supersymmetric 
particles. A large number of new parameters (sparticle masses, mixing angles, 
e.t.c.) appear in the calculations, enlarging therefore the parameter space and 
making difficult the consistency of the predicted branching ratios with the 
experimental bounds. However, in the context of unification and low energy 
phenomenology scenarios, these processes can provide useful constraints on 
the parameter space. 

We will briefly present the minimum number of inputs necessary to deter- 
mine all low energy parameters entering in a lepton flavour violating process. 
In the context of supersymmetric unified models, we assume a universality 
condition for the scalar masses at the unification scale My. The general for- 
mula, at this scale, is fhf(Mu) = (1 + qi)m\i 2 , where m 3 / 2 is the gravitino 
mass, % is a flavour index and qi is the modular weight of the corresponding 
field. This tree-level contribution is flavour diagonal. Non-diagonal terms are 
expected to appear through non-renormalizable terms with an expense of an 
extra parameter e, namely e = (0)/M where ((f)) is a singlet field vev and 
M a Planck-scale mass. The magnitude of the vev (0) can be fit from the 
fermion sector. A crucial role is also played by the gaugino soft masses, the 
trilinear soft parameter A as well as the Yukawa couplings X t and A;, at the 
unification scale Mu (we assume that Xb(Mu) = X T (Mu))- In the minimal 
scenario, the gaugino masses at the unification scale are determined in terms 
of the universal mass parameter mm. Thus, at Mu, we use a minimum set 
of parameters, namely (mu2, Tn.3/2, A*, A, X t , A&), together with the value of 
the common coupling olu and the unification scale Mu in such a way that 
after the renormalization group running we obtain a consistent set of all low 
energy measured quantities. For any acceptable such set, we calculate the 
branching ratios of the flavour violating processes. 

Figure 1 shows the one-loop diagrams relevant to the [i — > process. 
The corresponding r — > /i7-decay is represented by an analogous set of 
graphs, n — > 3e proceeds through the decay of the (now virtual) photon 
to an electron-positron pair. There are also box-diagrams contributing to 
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Figure 1: The generic Feynman diagrams for the /x — > e7 decay. / stands for 
charged slepton (a) or sneutrino (b), while and represent neutralinos 
and charginos respectively. 

this process, they are however relatively suppressed. 

The electromagnetic current operator between two lepton states k and lj 
is given in general by 

Tx = (h\( P - q)\Jx\lM) 

= Ufa- q){ mj ia x ^ (A L M P L + AfjP R ) + 

(g 2 7A - qxl • q) (A L E P L + A*P R ) } Uj (p) (13) 

where q is the photon momentum. The Am's and A^'s have contributions 
from neutralino-charged slepton (n) and chargino-sneutrino (c) exchange 

aL,R aL,R i aL,R aL,R aL,R. . /iL,i? (■\A\ 

M — A M(n) ~+~ ^Af(c)' 71 £ — ^(n) ^(c) l 1 ^ 

The amplitude of the process is then proportional to T x e x where e A is the pho- 
ton polarization vector. The easiest way to determine the loop momentum 
integral contribution to the A's is to search, in the corresponding diagram, 
for terms proportional to (p • e) and (q • e). The coefficient of the former 
is proportional to the momentum integral contribution to the <j\p term in 
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( |T3"D , while the coefficient of the latter is proportional to the difference of the 
momentum integral contribution between the <j\p and the (q 2, j\ — q\j ■ q) 
terms. Defining the ratio x = M 2 /m 2 , where M is the chargino (neutralino) 
mass and m the sneutrino (charged slepton) mass, the following functions 
appear in the A M term 

A-M{n) 



Am(c) 



where is the mass of the lj lepton, while for the Ae we have 

A E (n)- ( T J ^(2-9a; + 18a; 2 -18a; 3 + 6x 3 loga;) 

A E (c)- ^^(16 - 45z + 36x 2 -7a; 3 + 6(2 -3z) logs) ( l>) 

Notice the lack of terms proportional to the gaugino mass M which cancel. 
The Branching Ratio (BR) of the decay lj — > U + 7 is given by 

BR{h - kl) = ^ ((Am) 2 + (A*) 2 ) 

Our approach to determine the mixing effects is the following: 

• For each particular model we construct the lepton and left and right 
slepton mass matrices and determine the corresponding diagonalizing 
matrices. 

• From the input parameters at Mu (ajj, Mjj, mi/2, e.t.c.) and using 
the RGEs, we determine the soft masses for gauginos and sleptons, 
the Yukawa couplings At and A& and the /1 and A parameters. We are 
using two-loop-/? functions and incorporate threshold effects for the 
scalars and gauginos. In Table [3] we show the results of the RGEs 
running for four characteristic cases (In the runnings we have assumed 
the condition mi/2 > rn 3 / 2 . We have not explored the parameter space 
where the inverse relation holds). 

• Using the above values and the form of the mass matrices for gaugi- 
nos and sleptons we determine the corresponding mass eigenvalues and 

7 
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T (l — 6x + 3x + 2x — 6x log 2) and 



^(l-x 2 + 2xlogx)^- 

QMZ^i (2 + 3x — 6a; 2 + x 3 + 6x log x) and 
1 w (-3 + Ax - x 2 - 2 logx)^- 



(15) 



(l-x) 3 



eigenstates. The diagonalizing matrices determined in the first step are 
used to transform the slepton mass matrices in the desired basis where 
the lepton mass matrix is diagonal. 

• Having all relevant diagonalizing matrices (for charginos, neutralinos 
and sleptons), plus the eigenvalues of the gauginos and sleptons we can 
construct the amplitude of the process. 



Case 


(a) 


(b) 


(c) 


(d) 


INPUTS 


au\ Mu/10 w 


24.88, 1.0 


25.07,0.67 


25.00,0.67 


25.00,0.67 


m 3 /2, mi/2 


138, 350 


150, 380 


160, 390 


160,390 


At, At 


0.85, 0.10 


0.90, 0.10 


1.05, 0.08 


1.05, 0.05 


A, n 


-207, -395 


-225,-395 


-240,-395 


-240,-395 


RESULTS 


Mi, M 2 , M 3 


151, 295, 820 


166, 321, 865 


171, 329, 888 


170, 329, 890 


^(1,2) > ^(3) 


281, 279 


303, 301 


314, 313 


316, 314 


^(1,2)) ^(3) 


193, 187 


209, 203 


219, 215 


219, 218 


~ Q ~ Q 


750, 678 


789, 712 


811, 730 


813, 733 


m" 2 , rrig 


718, 564 


754, 589 


775, 598 


778, 599 


mf 2 , ^3 


718, 564 


754, 589 


775, 598 


778, 599 


mf 2 , ^3 


714, 199 


749, 219 


771, 244 


773, 252 


A T , fi 


10.3, -324 


8.3, -319 


3.3, -299 


6.9, -301 


tan/5, m t 


14, 185 


14, 188 


11,195 


7, 194 



Table 1: Inputs and outputs of the RGEs running for four representative 
cases (masses are in GeV). 



We present here the relevant mass matrices of two models whose successful 
fermion mass hierarchy is predicted by U(l) symmetries. As a first example, 
we use the scalar mass matrices obtained in a simple SU(3) x £77(2) x £7(1) 
model with an additional £7(1)/ symmetry ||. After the implementation of 
this symmetry, the fermion matrix for leptons in this model is given 

/ g2|a+6| ~\a\ ^\a+b\ \ 

m e w gl a l e 2 l b l el fe l m T (17) 
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where the parameter e is some power of the singlet vev scaled by the uni- 
fication mass, while a, b are certain combinations of the lepton and quark 
U(l) /-charges. Order one parameters in front of the various entries (not cal- 
culable in this simple model) are assumed to reproduce exactly the fermion 
mass relations after renormalization group running.^ The scalar mass matri- 
ces of this model are given in pj. For the sleptons we obtain 



(18) 



A successful lepton mass hierarchy in this case is obtained for the choice 
a = 3, b = 1 and e = 0.23. 

As a second example in the calculation of the rare processes, we apply 
our results in several textures obtained in context of the SU(4) model in 
ref.JTj. As far as the symmetric fermion mass textures are concerned, the 
analysis follows the same lines as in the first example, since the matrices 
obtained are similar to those above. The non-symmetric case is however 
completely different. There, the violations are much harder and the bounds 
on the various scalar mass parameters increase substantially. In order to be 
specific, we work out a particular example based on the charge assignment 
of Table |. The field assignment is F + F = (4, 2, 1) + (4, 1, 2) for the three 
generations, H + #(4,1,2) + (4,1,2) for the higgses and h = (1,2,2) for 
the bidoublet including the standard higgs doublets. Notice that this U(l) 



Ft 


F 2 


^3 


F x 


F 2 


^3 


h 


H 


H 


4 





-1 


3 


-1 


1 





1 


-1 



Table 2: U(l)f charges of fields in the model. 

symmetry is anomalous. However the mixed anomalies are zero. This fact 
allows for a Green-Schwarz anomaly cancellation mechanism. The fermion 
mass matrices are generated by operators of the form 

FFjh^e* 1 (19) 

where 5 = HH is an effective singlet generated by the higgs fields H, H. 
1 We use the Yukawa coefficients en = 4.0; 0x2 = c 2 i = 0.9; C22 = 1.08; C33 = 1.9. 
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Taking into account Clebsch-Gordan coefficients derived from these op- 
erators, as in ref. J7| and using the values e ~ 0.14, 5 ~ 0.21, we arrive at 
the following lepton Yukawa matrix at the unification scale 



/ v/27712 

V2r)2i 37722/ V$ 
V V33 



(20) 



The following unification scale input parameters give the correct lepton mass 
spectrum 



V22 



2.88 x 10~ 2 , 
V33 = 1-18, 



V12 

V32 



2.81 x 10~ 3 , 
7.28 x 10~ 2 . 



7721 = 1.30 x 10" 



(21) 



The scalar mass matrices are determined from the U(l) charges chosen 
to give correct predictions for the fermions. Since left and right fields have 
different charges, we obtain two different types of scalar mass matrices 



mf 



( 1 

e 4 

u 

/ 1 

e 4 

u 



e 4 
1 
e 

e 4 



e 5 



e 2 



(22) 



(23) 



The two expansion parameters e and e are defined as follows 



{6) /Mi 



u ~ e 



(9)/M b 



(24) 



In Table ^ we show the branching ratios for the four cases appearing in 
Table [I] and for the two models discussed above. We also give the L and R 
amplitudes for the neutralino and chargino exchanges for each case (see Eqs 

USE)- 

From Tables [L] and |3], we can conclude that lepton flavour violation in a 
viable class of supersymmetric unified theories puts non-trivial constraints on 
the scalar mass spectrum. In the case of symmetric textures of mass matrices, 
the n — > e7-branching ratio is found to exceed the present experimental 

l 3/2 



bound (4.9 x 10 n ) for values of the gravitino mass parameter m 3 /^ less than 
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CASE 


(a) 


(b) (c) 


(d) 


symmetric textures in U(l) /-models. 


A L 
aR 

A M{n) 

^Af(c) 
AR 

A M(c) 


3.58- 10- 11 
-9.47- 10~ 12 
4.13 • 10~ 14 
8.43 • 10~ 12 


2.74- 10- 11 
-6.79 • 10~ 12 
3.57- 10~ 14 
7.30 • 10~ 12 


1.22- 10- 11 
-2.04- 10~ 12 
1.44 • 10~ 14 
2.93 • 10~ 12 


1.83 • 10" 12 
-1.11 • 10~ 12 
8.84- 10~ 15 
1.79 ■ 10- 12 


BR(fj, -> e 7 ) 


1.03 • 10^ 1U 


6.01 • 10" 11 


1.20 • 10~ n 


3.07 -lO" 13 


M(n) 
M(n) 


9.15 • 10~ 10 
2.45 • 10~ 10 
-1.36 - 10- 11 
-2.26 • 10- 10 


7.05 • 10- 10 
1.83 • 10~ 10 
-1.18 • 10- 11 
-1.96 • 10- 10 


3.23 • 10~ 1U 
6.71 • 10~ n 
-4.74- 10~ 12 
-7.87- 10- 11 


4.98 • 10" 11 
3.18 • 10~ n 
-2.91 • 10~ 12 
-4.80- 10- 11 


Bi2(r - /i 7 ) 


6.50 • 10~ 8 


3.84- 10~ 8 


8.08 • 10~ 9 


1.96 • 10~ 1U 


Non-symmetric fermion mass textures 


A M{n) 
AR 

A L 

A M{c) 
AR 

A M{c) 


-4.55 • 10~ 10 
6.98- 10- 12 
-2.41 • 10" 14 
-4.87- 10~ 12 


-3.56 • 10- 10 
5.30- 10- 12 
-5.98 • 10" 14 
-1.21 • 10~ n 


-1.70 • 10~ 10 
1.72 • 10~ 12 
-1.22 • 10~ 14 
-2.43 • 10~ 12 


-2.55 • lO^ 11 
1.59 • 10~ 12 
-1.48 • 10^ 12 
-2.97- 10~ 12 


BR{n -»• e 7 ) 


1.66- 10" 8 


1.01 • 10" 9 


2.30 • 10" y 


5.21 • 10" 11 


%(n) 

4 R 

M(n) 
A L 
aR 

A M(c) 


-8.51 • 10~ 10 
1.59 • 10~ 10 
-1.53 - 10" 11 
-2.55 • 10~ 10 


-6.56 • 10- 10 
1.16 • 10~ 10 
-1.32 • 10~ n 
-2.21 • 10~ 10 


-3.03- 10- 10 
3.38- 10- 11 
-5.31 • 10~ 12 
-8.86- 10- 11 


-4.39 • 10" 11 
2.93- 10- 11 
-3.27- 10~ 12 
-5.40 • 10~ n 


BR(t -> /j 7 ) 


6.08- 10~ 8 


3.67- 10~ 8 


7.83 • 10- y 


2.27- 10~ 1U 



Table 3: Branching ratios for the two processes /i — > e 7 and r — > /i 7 , corre- 
sponding to the four cases of Table 1 and for the two models with symmetric 
and non-symmetric mass textures. We also show the L and R amplitudes 
with neutralino and chargino exchanges. 
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about 150Gev, in particular when tan (3 obtains intermediate or higher values. 
In our approach, the gaugino mass parameter 7711/2 is a ls° found to have a 
lower bound (> 350Gev) for consistency with the experimental /1 — > ey- 
bound. In the case of non- symmetric masses matrices, the bounds are even 
higher: the three cases (a, b, c) of Table 3 violate the present experimental 
bounds, implying therefore a rather heavy sypersymmetric mass spectrum. 
For small tan (3, however, we obtain results consistent with the experimental 
bounds as in case (d) . Actually, we see a strong dependence of the branching 
ratio on tan/3. As tan/? — > 1, the neutralino exchange processes get smaller 
and smaller. The chargino ones, while individually remain of the same order, 
they exhibit an increasing cancelation, rendering the branching ratio much 
lower that the experimental bound. The r — > /i7~rare decay (branching ratio 
< 4.2 x 10 -6 ), does not put further constraints as can been checked from the 
same Table. Moreover, the branching ratio for the /1 — > 3e-decay in all cases 
is found about three orders smaller than BR(n — > ej). 

Returning to the results of Table 1, we infer that non-observation of the 
— > e7-decay implies that all scalars appear with masses at least heavier 
than about 200Gev. Although our results are only for two specific examples, 
their main characteristics are rather generic and all sparticle mass bounds 
from flavour decays are expected always much larger than those obtained 
from other types of experiments. In this selected region of the parameter 
space (i.e. mi/ 2 > m 3 / 2 ) the sleptons are the lightest scalars. Clearly, the 
relatively large flavour violations are due the large tan /3-effects as well as to 
the fact that U(l) /-models imply also mixing effects to the scalar sector. The 
rather high scalar mass bounds could be considered as an ominous perspec- 
tive for models with large tan/3 and non-symmetric fermion mass textures, 
in particular for those who envisage a relatively light supersymmetric mass 
spectrum accessible to future experiments. We should note however, that 
slepton masses of this order are in the range of the LHC Indeed, slepton 
decays can in principle be detected in CMS with a mass up to 350Gev[^J 
Finally, we wish to comment that in more complicated structures with cyclic 
permutation symmetries between generations and universal anomalous U(l)- 
factors may prevent mixing effects in the supersymmetric mass matrices||. 
In such cases, the above constraints are relaxed. 

In conclusion, we have examined lepton flavour rare processes in a class of 
supersymmetric gauge theories where an additional U(l) symmetry discrimi- 
nates the three families. Such symmetries, are capable of generating success- 
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fully the hierarchical fermion mass spectrum and the Kobayashi-Maskawa 
mixing in the hadronic sector, when non-renormalizable contributions are 
taken into account in the superpotential. They, however, imply mixing effects 
through non-renormalizable terms in the Kahler potential, and consequently, 
in the scalar partners of quarks and leptons, leading thus to hard flavour vi- 
olations. These violations are even larger when models with non-symmetric 
textures and large tan/3 values are considered. As a result, stringent con- 
straints are found in the sparticle spectrum. Such bounds could be relaxed 
if, for example, additional symmetries of the Kahler potential are imposed, 
so that the scalar partners of fermions remain flavour diagonal, while at the 
same time they do not appear in the superpotential. We hope to come back 
to this issue in a future publication. 
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